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Abstract 



^ We elaborate how to construct the interweaving chiral spirals in (2+1) dimensions, 

'^ that is defined as a superposition of differently oriented chiral spirals. We divide 

I the two-dimensional Fermi sea into distinct wedges characterized by the opening 

CS| angle 20 and the depth Q ^ pf, where p^ is the Fermi momentum. Each wedge 

f""^ earns an energy gain by forming a single chiral spiral. The optimal values for O and 

O Q are chosen by the balance between this energy gain and the energy costs from 

, the deformed Fermi surface (dominant at large O) and patch-patch interactions 

:• (dominant at small O). We estimate these energy gains and costs by means of the 

• i-H expansions in terms of l/A/'c, Aqcd/Q? and O using a non-local four- Fermi interac- 

/\ tion model: At small l/N^ the mass gap (chiral condensate) is large enough and the 

^ interaction among quarks and the condensate is local in momentum space thanks 

to the form factor in our non-local model. The fact that patch-patch interactions lie 

only near the patch boundaries makes it simple to embed many chiral spirals. We 

identify the dominant and subdominant terms at high density and formulate the 

perturbation theory in terms of subdominant terms suppressed by Aqcd/Q ot O. 

The kinetic term in the transverse directions turns out to be subdominant, so that 

we can utilize techniques in the (l+l)-dimensional system. Our central results at 

the leading order of the 1/A^c and the Aqcd/Q expansions are that the total energy 

gain is ^ pfAqqj^ with O ^ (Aqcd/pf)^^^- Since O decreases with increasing pF, 

there should be phase transitions associated with the change in the wedge number. 

We also argue the effects of subdominant terms at lower density where the large- A/'c 

approximation is more reliable. 
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1 Introduction and Our Central Results 



In recent works [1,2,3,4,5,6], it has been argued that there is a new state of 
QCD matter at high baryon density and low to intermediate temperatures. 

This novel state is called Quarkyonic matter since it exists at densities very 
large compared to the QCD scale so that the Fermi sea is best thought of in 
terms of quark degrees of freedom, but nevertheless is confining. It may be 
thought of as a Fermi sea of approximately free quarks, but with thermal and 
Fermi surface excitations made of color-confined mesons and baryons. The 
name "Quarkyonic" expresses this Yin- Yang dualism. 

While the arguments for the existence of this matter are rigorous only for a 
very large number of colors, this may not be such a bad approximation for 
some range of density for A^c = 3. The inter-quark potential inferred from the 
charmonium spectrum is linear out to distances of ^ fm at least indicating 
that the quark-pair production is not very efficient in tempering its growth. 
One way to understand this is the large- A^c limit where quark pairs are the 
source of 1/A^c corrections [7]. Also in lattice Monte-Carlo studies, the (pseudo- 
)critical temperature of the QCD phase transitions is a slowly varying function 
of baryon density, as long as densities are sufficiently small to estimate such 
effects reliably [8]. 

In this environment of high baryon density, one might expect that chiral sym- 
metry would be restored while the nature of quark confinement survives. In 
fact for a spatially homogeneous chiral condensate, several computations have 
confirmed this expectation [3,4]. Such characterization is, however, challenged 
by later analysis [9] and, moreover, simple phenomenological arguments sug- 
gest that chiral symmetry be broken in a confining phase of QCD [10]. 

Chiral symmetry restoration for a homogeneous phase is understood as follows. 
In a homogenous phase, the scalar mesons that condense to form the chiral 
condensate have zero net momentum. A usual chiral condensate composed of 
quarks and anti-quarks is not energetically favored since to create anti-quarks 
costs a large energy of the order of //q c^ p^ (where /Xq and p^ are quark 
chemical potential and the Fermi momentum, respectively) when they are 
popped up from the Dirac sea to above the Fermi sea. An alternative way to 
have a homogeneous chiral condensate is to pair up quarks with quark-holes 
near the Fermi surface; see the left panel in Fig. 1. In the presence of the Fermi 
sea, quarks must have finite momenta, and holes also must have momenta of 
the same order but with an opposite sign to make the net momentum of scalar 
mesons zero. Therefore, the relative momentum of the quark and the quark- 
hole in a pair is enormous, and in a confining theory, the string tension would 
require that the excitation energy of such a bound state be of the order of 
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Fig. 1. The particle- hole pairing with the confining interaction. (Left): The homo- 
geneous chiral condensate with total momentum ptotai — 0. The relative momentum 
between a particle and a hole is large. (Right): The inhomogeneous chiral conden- 
sate with ptotai — i2pF^z- The relative momentum is small. A superposition of 
pairs with momenta ptotai — '^Pf^z 9^nd ^ptotai — —'^Pf'^z creates chiral spirals. 

2/Xq relative to that of the scalar meson in vacuuiiJM It is unlikely that the 
condensation of such highly excited scalar mesons earns the energy gain, and 
so chiral symmetry restoration is concluded. 

It is, however, possible to have particle-hole condensation in the form of charge 
density waves [11,12,13,14,15], which is a similar state to the p-wave pion 
condensation in nuclear matter [16]. Early studies using perturbative gluon 
propagators [12,13] argued that the charge density waves do not realize unless 
the number of colors is huge. These arguments, however, did not take account 
of the confining effects. More precisely speaking, the attractive force in the 
infrared (IR) sector was not strong enough to overtake a screening effect. In a 
recent paper [17], several of us have argued that in Quarkyonic matter transla- 
tionally non-invariant chiral condensates form as chiral spirals. The argument 
for a translationally non-invariant condensate follows again from a particle- 
hole pair near the Fermi surface; see the right panel in Fig. 1. A difference 
from the homogeneous condensate is that quarks and quark-holes co-move in 
the same spatial direction, and they exchange only small momenta of the order 
of Aqcd, where Aqcd denotes the typical scale of QCD. Therefore forming a 



^ Presumably this is a sufficient condition not to have the homogeneous particle- 
hole condensation. Actually, even without the confining force, it is likely that the 
density- wave type is more favored. See discussions below Eq. (43). 
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Fig. 2. The two-dimensional slice of the Fermi sea, suggested in Ref. [19]. Rotational 
symmetry is spontaneously broken from continuous to discrete one. The number of 
patches accompanying chiral spirals in different directions increases with increasing 
density. 

bound state does not cost much energy in contrast to the homogeneous case, 
so the condensation is possible. One finds that the optimal mode of conden- 
sation is a linear combination of the chiral condensate, ('^'0), and a spin-one 
isosinglet odd-parity excitation, {ipa^^ip)^ where z is the direction of motion 
of the wave [12,17], which will be referred to as the "longitudinal" direction 
and the orthogonal are called the "transverse" directions throughout this pa- 
per. The chiral spiral is characterized by the spatial oscillation between these 
two modes. This combination can be naturally interpreted as a superposition 
of particle-hole condensates with momenta ^ ±2pYnz- At high density, like 
heavy-quark symmetry, there emerges approximate SU{2Nf)^ x SU{2Nf)_ 
symmetry [12], where ± expresses (l+l)-dimensional chirality that charac- 
terizes the moving directions along the z axis. Therefore, after the formation 
of chiral spirals, there are (2A^f)^ Nambu-Goldstone (NG) modes: (2A^f)^ — 1 
isospin-spin excitation modes and one phonon mode associated with sponta- 
neous breaking of spin-chiral and translational symmetry, respectively^] These 
results were derived by the dimensional reduction from the (3+l)-dimensional 
self-consistent equations to those in the (l+l)-dimensional 't Hooft model for 
degrees of freedom near the Fermi surface. 

The stability analysis of Ref. [17] showed that Quarkyonic matter in the ab- 
sence of a chiral condensate was unstable with respect to the formation of 
a (l+l)-dimensional chiral spiral. Furthermore, it was suggested that many 
chiral spirals of different spatial orientations interweave to form a more com- 
plicated condensate [19]. This corresponds to a transition from a spherical 
Fermi surface into patches, inducing breaking of continuous rotational sym- 
metry down to discrete one (see Fig. 2). The number of patches increases as 
the density grows up and such phase transitions continue to occur until the 



^ The formation of a single chiral spiral breaks rotational symmetry in addition 
to translational symmetry. Since the translation and rotation are not independent, 
only one phonon mode appears as an NG mode [18]. 
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Fig. 3. The three dimensional plot of the phase diagram in jiq-T-E space {£ is the 
energy density of the system). Chiral symmetry in Quarkyonic matter is broken by 
the formation of the interweaving chiral spirals (ICS). The stair-like growth of the 
energy density along fiq axis expresses the discontinuos change in the shape of the 
Fermi sea (Fig. 2). 

screening effect on gluons strongly reduces the IR attraction between a pair of 
a quark and a quark-hole. Such reduction happens around the density scale, 

/i,^iVV2AQCD[l]. 

These behavior are sketched in a the phase diagram in fiq-T-E space, as shown 
in Fig. 3 {8 is the energy density of the system)]^ Quarkyonic matter starts 
to appear just above jiq ^ Aqcd, where a transition from nuclear to quark 
matter quickly occurs, and continues to exist up to jiq ~ A^^^/^Aqcd- In the 
Quarkyonic region, the stair-like growth of the energy density along /x^ axis 
reflects the discontinuos changes in the shape of the Fermi sea (Fig. 2). At 
larger /i^, the screening of the IR attraction reduces the size of the chiral spiral 
condensate, and accordingly, the interval of stair in jiq axis and jumps in £ 
become smaller. The shape of the Fermi sea smoothly approaches spherical 
one. 

While Refs. [17,19] have argued for chiral spirals using the conflning interac- 
tions, several aspects have not been explicitly shown due to technical difli- 
culties and/or conceptual uncertainties in treating the deep IR structure of 
confining forces. On the other hand, while we postulate that the confining 
force could give a sufficient condition to drive chiral symmetry breaking, it is 
certainly not a necessary condition. The aforementioned phenomena may ap- 
pear in a wider class of models that encompass chiral symmetry breaking even 



^ Here we consider chiral limit for the light flavors and ignore the electroweak 
interactions. 
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Fig. 4. The leading self-energy diagram at zero density. At large A/'c, we have only to 
keep the rainbow ladder. (Left) The diagram in terms of QCD dynamics. Integrating 
the temporal component out, we can interpret the loop with momentum k as the 
condensate, {ilj(k)ilj(k)), which is made of particle-antiparticle with momentum k. 
(Right) The corresponding diagram in our model. The soft-gluon exchange part in 
QCD is replaced with the form- factor function whose strength damps as k and pgo 
far apart. The diagrams for the 1/Nc corrections will be given in Fig. 17 in Sec. 8. 

without confinement. Indeed, what is relevant for interweaving chiral spirals 
are Fermi surface effects and the IR enhancement of the interaction, but not 
precise knowledge about the deep IR region. Taking this viewpoint, we will 
characterize chiral symmetry breaking at high density by a simple, tractable 
model in which we can explore analytic insights. 

We will use an effective field theory to describe QCD at large Nc. This model 
will be apparently similar to the Nambu-Jona-Lasinio (NJL) model [20,21] in 
which the dominant interacting process in the large- A'c limit is the scattering 
of particles and the condensate. The crucial difference of our model from the 
usual NJL model is that the interaction vertex has a form factor that mimics 
the IR enhancement of the non-perturbative gluon propagator. We denote a 
form- factor scale of the model as Af (^ Aqcd), beyond which interactions are 
negligible compared to the IR interaction. 

Form- factor effects lead to the following consequences at large Ac- The in- 
teraction between a quark and a condensate becomes strongly momentum 
dependent; see Fig. 4. They decouple one anotheip] if the relative momentum 
between the quark and the condensate is much larger than Af , reflecting com- 
posite nature of the quark condensate. As a consequence, the quark mass gap 
damps when the quark momentum is far away from the domain of condensa- 
tion. In vacuum, the damping scale of the mass function Ac may play a similar 
role to the ultraviolet (UV) cutoff, Anjl, in the usual NJL model (see the left 
panel in Fig. 5). In this sense, the form factor can naturally remove the UV 
cutoff artifact of the NJL-type model. 



^ This sort of picture has been discussed for the high-lying mesons and baryons 
[22]. See also Ref. [23] for some caveats on this picture. 
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Fig. 5. Schematic figures for the mass gap function. (Left) The zero density case. 
(Right) The finite-density case with chiral symmetry breaking near the Fermi sur- 
face. 

All of these aspects are crucial when we consider dense quark matter. Since 
condensation phenomena should happen near the Fermi surface, the effective 
UV cutoff should appear in the distance from the Fermi surface, not that from 
vacuum (see the right panel in Fig. 5). If we fixed the effective UV cutoff to 
be the same as the vacuum value by hand, arguments on chiral symmetry 
breaking would not make sense for /Xq > Anjl- Indeed, because of Pauli block- 
ing, phase space for quarks contributing to the condensate would disappear, 
leading to chiral symmetry restoration as a cutoff artifact. Therefore it would 
be desirable to derive an effective UV cutoff dynamically for each density. We 
claim that the introduction of the form factor gives a natural extension of the 
treatment of the zero-density NJL model to that at finite density. Also we 
make a remark that the so-called Debye cutoff frequency is introduced in this 
way from the Fermi surface in the standard BCS theory. 

With this modified NJL model at hand, the purpose of this paper is to give 
detailed and analytic insights on the interweaving of chiral spirals and on the 
associated breakdown of the rotational invariance. A key observation in our 
model treatment with a form factor will be that at large Nc^ the particle- 
condensate interaction (i.e. the interaction of particles scattering off the con- 
densate) happens locally in momentum spacq^ which will allow the system 
to simultaneously embed many chiral spirals at sufficiently high density. 

For the sake of simplicity, we will work in (2+1) dimensions, where the original 
Fermi surface is a circle and takes a simple geometric structure even after for- 
mation of many chiral spirals. An extension of this study to higher dimensions 
might be technically difficult but conceptually straightforward. 

We approximate the patched Fermi surface by a polygon of degree A^p. In Fig. 6 
we show the two-dimensional Fermi surface and some polygon approximations 
to it. We will look for an energy minimum at a non-zero value of Ap, and we 



^ The 1/Ac corrections (as shown in Fig. 17) will violate this locality; see discus- 
sions in Sec. 8. 






Fig. 6. The leftmost is the Fermi circle in spatial two dimensions. The middle fig- 
ure is a square approximation to the circle. The wedges are shown corresponding 
to discrete sub-area of the circle. The rightmost figure is a higher-order polygon 
approximation to the circle. 

will find that there exists such a minimum that depends upon density. We 
can think of each sub-sector that constitutes the polygon shape as a wedge. 
The wedge is characterized by an opening angle 20 and a depth Q as shown 
in Fig. 7. The depth Q will be of the order of the Fermi momentum pp, and 
the opening angle of the wedge is constrainecp] by 2 x 2A^p0 = 27r. The 
surface thickness Apermi ("^ Af) characterizes the momentum scale for which 
non-perturbative Fermi surface effects are important. 

The use of the wedge shape is motivated by the following reasons: In order 
to maximize the energy gain from condensation effects, each patch should 
contain only one chiral spiral by aligning total momenta of a bunch of particle- 
hole pairs. If we had a misalignment, interplay among chiral spirals would 
reduce the size of the gap, as exemplified in Sec. 7.3. So we have to look 
for the most effective shape to achieve the alignment of total momenta. An 
obvious candidate is the flat Fermi surface with which particles and holes 
participating in the condensate can stay close to the Fermi surface, saving 
the virtual excitation energies. Other shapes require some of particles and 
holes with larger excitation energies, so are not effective to create a bigger 
condensate. 

What is the principle to determine the number of patches? It is essentially de- 
termined by the balance between the kinetic-energy cost and the condensation- 
energy gain. Without condensation effects, the Fermi surface would simply 
take a circle shape that minimizes the kinetic energy. Once condensation ef- 
fects are turned on, the associated energy gain can overtake the kinetic-energy 
cost, changing the shape of the Fermi sea from a circle to a polygon with dis- 
crete number of wedges. The condensation effects appear in two places. One is 
the energy reduction of single-particle contributions inside of one wedge, as is 
the case in the single chiral spiral problem. The other is the inter-patch inter- 



^ Here a factor 2 is included since we will take one patch as a set of one wedge 
and the other wedge in the opposite side of the Fermi sea. See Sec. 3 for details. 
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Fig. 7. The basic wedge from which the polygon is constructed. The wedge has an 
opening angle 20, a depth Q, and a surface thickness Apermi- 

actions among different chiral spirals. The latter will provide an energy cost, 
that is, chiral spirals destroy one another if their wavevectors are different (as 
explained in Sec. 7, and also some results from the existing literatures will be 
addressed in Sec. 8.1). This means that we should divide the Fermi sea into 
not too many wedges to get the largest energy gain. Therefore the size of one 
wedge tends to be as large as possible until the kinetic-energy cost becomes 
too big. 

Let us outline how to optimize the shape of the Fermi sea. In the following, 
we will consider the canonical ensemble, in which the particle number is fixed. 
Then the Fermi volume must be conserved, so it follows that 



j9^e = g2tane. 



(1) 



Here the LHS is a fermion number for the plain circle shape of the Fermi sea, 
while the RHS is for the deformed Fermi sea. For small B we can approximate 
the RHS by the Taylor expansion in terms of B, with which we can solve Q{0) 



as 



Thus, the size of one patch in the transverse direction is 



(2) 



(one patch size) = Q(B) tan B = p^Q (l + -B^ + — B^ + •••)• (3) 

^ 6 360 ^ 

Now that Q is solved as a function of B, the multiple chiral spiral states are 
characterized by pp, B, and the single-particle mass gap, M. We will perform 
the energy minimization by taking B and M as variational parameters. 



Provided that the energy minimum exists for B ^ 1, we can expand the 



energy density by powers of 9. It will turn out that the expression takes the 
following form, 

£-i{M) 



e 



+ So{M) + 82{M)Q' + f4(M)e^ + . . . . (4) 



The 1/B term can appear simply because A^p = 27r/B, but it will disappear 
for vanishing M, that is, ^_i(M ^ 0) = 0. Expanding Sn{M) by powers of 
M/py^ we expect 

S^{M) = c^^^pI + c^^^pIM + cf ppM^ + • • • , (5) 

where we did not write possible non-analytic terms explicitly. (In what follows, 
solving the gap equation, in fact, we can arrive at the above form of the 
expression.) At sufficiently high density, one has only to keep the leading term 
in the M/pp-expansion for each Sn{M). 

Let us first discuss terms insensitive to details of condensation effects, comput- 
ing at M = 0. The first non- vanishing contribution of 0{py) with non-trivial 
G dependence should arise from the kinetic-energy cost for the deformation 
of the Fermi surface. The energy at M = is an increasing function of B, 

where the first term gives the trivial contribution which should be subtracted 
out. It is extremely important to notice that the non-trivial deformation energy 
does not appear until O(B^). The volume conservation of the Fermi sea cancels 
the B^-term out from the average kinetic energy. 

Terms beyond (9(B^) are much smaller than p^B^ and irrelevant in our mini- 
mization procedure. Thus, the energy minimum will be found by balancing the 
PyQ"^ term with condensation effects which yield terms with smaller powers of 
B than the p^Q^ term. In the following we concentrate on the estimation of 
such condensation terms. 

The condensation effects depend on interaction properties of models. The point 
of our model is that at large Nc the single-particle dispersion of a fermion with 
momentum p is determined only by condensates within the domain of the size 
^ Af around p. It means that if we consider the gap for single particles farer 
from patch boundaries than Af, it will be affected only by the single chiral 
spiral, not by adjacent chiral spirals. 

The above argument suggests that we have to evaluate the mass gap differently 



10 



depending on the different domains of 9. then there are intersection points 
of more than one patch that have interactions among them. To solve the 
gap equation, hence, we have to take into account the influence of several 
chiral spirals simultaneously. This is a rather technically complicated problem. 
Fortunately, the energy minimum in our problem will be outside of this G 
domain. 

We can self-consistently show that the transverse size of one patch is much 
larger than Af, i.e. 

Pf© > Af . (7) 

Once this condition is satisfled, the single-particle gap in one patch can be 
determined independently from B, except in the region close to the patch 
boundaries. We denote such a solution as M = Mq ^ Af. Then the energy 
gain from condensation effects should be 

(energy gain) - A^p • A^c (Af Pf tan 9) Mq - N,Mo Af pp [ 1 + y + • • • ) , (8) 

where Af pp tan B is the one-patch phase space where the condensation occurs. 
One important observation here is that, while the gap is insensitive to B, the 
phase space has B dependence, so that the leading term is B independent 
after multiplying a patch number Ap. 

Let us see contributions at the intersection point of two adjacent patches. The 
point is that a particle from one patch and condensates from other patches 
interact within a limited domain of ~ Af near the intersection points. Its phase 
space is independent of B. Therefore the contribution from the intersection 
points is 

(energy cost) - N^ ■ TV^A^ /(Mb) ^ ^ • Af /(^b) , (9) 

where /(Mb) is some function of the order Mb ^ Af with Mb be the mass gap 
near the boundary, and vanishes as Mb -^ 0. The contribution must be an en- 
ergy cost. The reason is that in Eq. (8) we overestimated the energy gain which 
should be reduced around the patch boundaries. The misalignment of chiral 
spiral wavevectors tends to destroy the different chiral spirals one another, 
and reduces the size of the gap at the intersection points. Diagrammatically, 
this contribution will appear as interactions among chiral spiral mean flelds 
in different domains. The presence of this term becomes more important for 
smaller B. 

Now we can express the energy density as a function of M and B. In the 
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Fig. 8. The condensation region near the Fermi surface. The thickness of the region 
in the radial direction is ^ Af. In the boundary region with the transverse size of 
^ Af, inter-patch interactions between the nearest-neighbor chiral spirals destroy 
condensates one another, reducing the energy gain from condensation effects. 

domain Af/pp <C B <C 1, it reads 

SS{M,e) 

- ^c(^^^§^ - coMoAfPF - csMoAfppe^ + C4j9^e^ + • • • y (10) 

where coefficients Cq, C2, • • • are positive, and we have subtracted the free Fermi 
gas contribution. The energy balance is schematically illustrated in Fig. 9. To 
assure this expression by microscopic calculations is our goal in later sections. 

When Af/pF < (Af/pp)"^^^ < B <C 1 is satisfied, the ppB^ term dominates over 
other terms, and the deformation energy supersedes the condensation energy. 
Thus, there is an upper bound of G, and the energy minimum should lie in 
the region, Af/pp <C 6 <C (Af/pp)"^^^. On the other hand, the lower bound 
of B will be set by the patch-patch interactions proportional to 1/0. In the 
region of current concern, we have 



d5S{M,e) 



- Ncl-^+^c^p^vOl] ^ 0, (11) 



e=eo V ^0 



for the energy minimum neglecting other terms. Therefore we find 




Go - - . (12) 



As we promised, we can confirm in this way that ppBo ^ (]9p/Af)^/^Af ^ Af, 
which surely justifies Eq. (7) posteriori 
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Fig. 9. The schematic energy landscape as a function of O. The region O <C Af/pp 
is beyond the apphcabihty of our analysis. 

The total energy is dominated by the condensation term ^ — ppAf , which is 
independent of Bq. The leading corrections come from the l/Oo term (patch- 
patch interactions) and the Gq term (the deformation energy) which are sup- 
pressed by (Af/pp)^^^ compared to the leading contribution. This implies that 
after solving the one-patch problem at the mean-field level, other patch-patch 
contributions can be treated within perturbation theory. We will formulate 
our perturbation theory within the domain of (Af/pp) <C B <C (Af/pp)"^^^- 

This paper is organized as follows. In Sec. 2, we define our model and describe 
consequences of form-factor effects in vacuum. In Sec. 3, the Fermi sea is 
decomposed into several segments. We formally separate the Lagrangian into 
the one-patch and the patch-patch interactions. In Sec. 4, the mean-field treat- 
ment for the one-patch problem is discussed. We first identify the dominant 
terms within one patch, and construct the mean field for chiral spirals as well 
as the mean-field propagators for quasi-particles. In Sec. 5, we give a formal 
expression of the perturbative expansion. In Sec. 6, we treat corrections from 
subdominant terms in one patch. It will be shown that subdominant terms 
are suppressed by powers of B or A^/Q. In Sec. 7, the inter-patch interactions 
at the patch boundaries are discussed. The size and sign of 1/B terms are 
estimated in both perturbative and non-perturbative manners. In Sec. 8, we 
argue possible impacts of several corrections ignored in this paper. We also 
review other works and attempt to place this work in perspective. Section 9 
is devoted to a summary and possible future directions. 



2 A Model: The Four-Fermi Interaction with Form-factor Effects 



In this section, we introduce explicit form factors into the NJL model, so 
that the interaction is cutoff when the momentum transfer becomes large. 
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This makes the NJL model non-local. (For other attempts to introduce non- 
locality which emphasize different aspects from ours, see Ref. [24].) We will 
see that effects of such a cutoff determine the momentum domains where 
condensation phenomena are relevant. This aspect is particularly important 
when we consider the very large Fermi sea in which condensation phenomena 
occur as Fermi surface effects, rather than the vacuum effects. 

Let us consider the scalar-scalar type of the four- Fermi interaction, 

fd?x{^^{x)f= fdx'f {i>{p + q)Hp)){i>{k)ijik + q)), (13) 

where we define a shorthand notation, 

r ^ r dg dpdfc 

W"^ (27r)6 ' ^^^^ 

and we did not explicitly write the coupling constant and Xq dependence of 
fermion fields. Since the four-Fermi interaction is not renormalizable beyond 
(1+1) dimensions, we need introduce some UV cutoff. 

We regularize the UV interaction by including form- factor effects, 

f d'x [i>ij{x)y ^ /dx°/ {i'{p + q)i^{p}){i'{kmk + q))ep,k, (15) 

where 

e^^, = e(K}-{p-kf). (16) 

This mimics the form-factor effects in large- A^c QCEQ and removes the UV 
divergences associated with interacting processes. 

The large- A^c QCD is mimicked as follows. The one-gluon exchange including 
non-perturbative effects are shown in Fig. 10(a). Its strength damps as the 
momentum transfer becomes large. We roughly take into account this prop- 
erty by introducing a step function, ^(A^ — {p — fc)^), keeping the interaction 
strength constant. In QCD, the cutoff scale Af should be taken to be of the 
order of Aqcd- 



^ Here we put the cutoff on the spatial- momentum p^, not on the Euclidean mo- 
mentum p\. So, results in this work are connected to those of large- A'c QCD in the 
Coulomb gauge, in which the dominant non-perturbative part is of the instantaneous 
type [25] . An alternative choice would be to introduce a cutoff on p\ keeping mani- 
fest Lorentz invariance. Such a treatment should mimic, for instance. Landau-gauge 
results in Euclidean space. However we do not know their Minkowskian behavior, 
so we have to compute quantities in Euclidean space. Then the price we have to 
pay is that a simple physical intuition does not necessarily work, especially at finite 
density. 
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Fig. 10. (a) The non-perturbative gluon exchange which is supposed to damp 
quickly in the UV region, (b) The color line representation of the one-gluon ex- 
change, (c) Our effective four- Fermi interaction including form factor effects. 

In Fig. 10(b), we show the color line representation to illustrate how the 
one-gluon exchange interaction should be contracted into a four-Fermi type 
interaction. Taking into account features in Figs. 10(a) and 10(b), we arrive 
at a simple model described in Eq. (15) and Fig. 10(c). 

A relevant consequence of our choice of the form factor is that the coupling 
between fermion fields and the condensate becomes strongly momentum de- 
pendent. In particular, fermions decouple from the condensation effects if they 
belong to domains with a momentum difference of the order of Af . 

Let us first see this property in the vacuum case by investigating the Schwinger- 
Dyson equation for quasi-particles. After picking up the pole, we havq 

where T>jn{k) is the self-energy and e{k) = yk'^ + S^(fc) is the quasi-particle 
energy. For its diagrammatic expression, see Fig. 4P1 When \p\ is very large, 
\k\ must be as large because of the form factor. In the asymptotic region 
|p] :» S(p), the Schwinger-Dyson equation looks like 

1 r dk ^ 

from which one can show that T>rn{p) damps faster than l/e{p) ^ l/|p] in 
the UV region, meaning that the scattering between quasi-particles and the 
mean- field condensates diminishes for large momentum. 



^ We ignored the Lorentz vector self-energy for the sake of simplicity. This simpli- 
fication should not alter our basic statements in this paper. 

^ The Fock-type gluon exchange, which is dominant in large- A'c QCD, corresponds 
to the Hartree term of the four-Fermi interaction model, so we can apply most of 
techniques used in the NJL model calculations. 
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This UV behavior tremendously simphfies considerations on the energy benefit 
from condensation effects. Due to decouphng from condensation effects in the 
UV region, the normahzed energy density, in which the energy density without 
condensation are subtracted out, is dominated by the IR contributions. Also, 
the couplings between the IR and the UV regions are allowed only in the 
limited region constrained by the form factor effects. Thus, we can proceed 
with our calculations independently from details of the UV physics. 

At finite density, the same argument applies if we replace the Dirac sea with 
the Fermi sea, namely, the excitation energy e{p) is measured from the Fermi 
surface instead of the vacuum. Then it follows that dominant contributions 
to the condensation come from fermions near the Fermi surface rather than 
those near vacuum. 



2. 1 Bosonization of the Model 



Following usual treatments, let us introduce the auxiliary fields and formally 
rewrite the action. Although this formal treatment is not mandatory, it would 
be practically convenient to use several relations derived in this framework. 

Since our four- Fermi interaction is not of the separable type, we have to modify 
conventional treatments slightly as 

/d^Y^ ^{Q)^{-q) -^Jdx'l^^^^q;p} 6,^, $(-<f;fc). (19) 

Here we did not write x^ dependence of the auxiliary boson field, $, explicitly 
for notational simplicity. Next we shift the boson fields to generate a four- 
Fermi interaction, 

$(g;p) ^ $(g;p) - i^{p + q)^{p) , (20) 

then, by adding these trivial terms in $, we can eliminate the four- Fermi 
interactions out from the original action to obtain the Yukawa-type model. 
The form of the vertex is 

/dx"/ ^{q-p)ep,ki'ik)ij{k + q). (21) 

We can obtain the gap equation from = /P$P^P^ ^^^^^e^*^ and find. 



/ 



^ 0,,,{m-q;k))-mmk + Q))} = O. (22) 



Since the equation should hold at arbitrary p, we conclude that ($(— g; k)) = 
{ip{k)ilj{k + q))^ which is independent of x°. 
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The relation between the mass self-energy and $ can be read off from the 
coefficient of ipip in Eq. (21), i.e. 



S™(g; k) = Jj^ Hq; P) e,,k . (23) 



In the stationary phase approximation of the ^-integral (which is rigorous in 
the large- A^c limit), we can replace $ with ($). One can check that substituting 
the gap equation for $ and taking g = (for uniform condensates) reproduce 
the Schwinger-Dyson equation in Eq. (17). 

Instead of computing self-consistent solutions precisely, let us consider what 
kind of Ansatz would capture basic properties of stationary solutions of $((f; p) . 
We will apply such arguments to more complicated case at finite density where 
condensation effects are important near the Fermi surface. 

In the vacuum case, we already know qualitative behaviors of S^((f;p) at high 
momenta from the previous arguments. Its value damps as |pj ^ oc, leading to 
damping of $(g;p), as seen in the gap equation. If we denote the characteristic 
scale of damping as Ac, presumably the simplest Ansatz would be 



<^{q-,p) = ^-d{Kl-p')5{q). (24) 

See also Fig. 5. In principle, not only A, but also Ac should be treated as 
variational parameters as functions of Af that is an intrinsic scale in our model. 

Perhaps it might be helpful to mention the relationship between this Ansatz 
and usual spatial-momentum cutoff scheme. In our model, conventional treat- 
ment of the NJL model is recovered by taking Af ^ cx) with Ac kept fixed 
to be the NJL-model cutoff Anjl ^ 600 MeV which semi-quantitatively de- 
scribes low-energy properties of the hadron phenomenology. Such successes of 
the NJL model imply that once we correctly cutoff the domain of condensa- 
tion^ we should be able to express the low-energy dynamics quantitatively. In 
our model, this cutoff scale should be dynamically determined by adjusting a 
value of Af . 

Adopting this interpretation, the effective cutoff of the NJL model should be 
measured from the Fermi surface, not from zero momentum, because conden- 
sation effects mainly appear near the Fermi surface. This leads to a picture 
such that chiral symmetry is restored deeply inside of the Fermi sea but is 
broken near the Fermi surface. 
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3 Decomposition of the Lagrangian into Multiple Patch Domains 



In this section, we decompose the NJL Lagrangian into different segments. 
For the sake of simphcity, we begin with the Lagrangian with one flavor, 



£ = ^ i # + ^ ((^^)2 + (^ i^^^)2) ^ (25) 

where we exphcitly factor out the A^c dependence of the interaction, so that 
G = 0{N^). In the (2+l)-dimensional system, G is dimensionful, and we take 
it to be ^ A^-^. Our Lagrangian has continuous chiral symmetry, U{l)i, x 
t/(l)R. In this work we ignore the U{l)x problem that is of 0(1/A^c) and thus 
neghgible. For the moment, we wiU not write the form factor exphcitly for 
notational simplicity. 

We introduce the unit vectors n^y which point to the center of the i-th patch, 
and the unit vector n^^ which is orthogonal to n^y. We project out the spatial 
components of vectors generically as 

Pi\\ =ni'p, pi± = ni± • p. (26) 

With this deflnition, fermion flelds can be decomposed into A^p momentum 
domains. 



2 P poo roo rpi^^tdiuG 

= T^T^ E / dpo / dpm / dpi± i){p) e-'P-^ 

(^ZTTj ^ J —oo J —oo -'— p^ijtanB 

^T.M^)- (27) 

i 

It is worth mentioning here that our formulation has potential relevance to 
analytic approaches based on the high-density effective theory [26]. A single 
patch includes a set of two domains with p^y > and p^y < for the positive 
energy states, in a similar way as the (l+l)-dimensional problem. Thus, the 
angle for one patch is 2 x 2B, and B should satisfy 

49 A^p = 27r . (28) 

Let us flrst decompose the free part of the Lagrangian. Since it is diagonal in 
momentum space, we have 

£'^- = ^t/;,^^Vi = E4'^ (29) 
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The decomposition of four- Fermi interactions O^^ are much more cumbersome, 
since the interaction terms combine different domains. Exphcitly writing, it 
foUows 

It can be decomposed into four types of interaction terms, 

i ^ i 

^ £g, = aT E ((V^i^i)(^i^/c) + (t/5a75V'i)(V^ii75^* 

Except the first hne, interaction terms involve fermions belonging to different 
patches. Now we have the Lagrangian separated into one-patch and patch- 
patch interactions, C = Z^i A"^^*^ + ^>C^^S where 

The one-patch Lagrangian will play a dominant role for condensation effects. 
We will first solve the one-patch problem, and then include contributions from 
other patches perturbatively. 



4 The One-Patch Problem at the Mean Field Level 



In this section, we provide a formal treatment of the one-patch Lagrangian. 
We consider sufficiently high density for which the interaction scale ^ Af is 
much larger than the transverse kinetic energy near the Fermi surface which 
is suppressed as ^ P±/Q where Q ^ Pf- Then the fermion dispersion relation 
is robust and does not rely on the transverse momentum; corrections from the 
higher spatial dimensions than the (l+l)-dimensional problem are generally 
suppressed by extra powers of l/Q. Therefore, the gap equation for the mean 
field effectively becomes (1+1) dimensional. We here analyze this quasi (1+1)- 
dimensional problem in detail. 
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We first summarize some convenient notations used in (1+1) dimensions. They 
are useful to identify the dominant and subdominant terms for the chiral spiral 
formation. Then we bosonize the dominant part of the four- Fermi interactions, 
and construct the mean-field Lagrangian for the chiral spirals. We prepare the 
mean-field propagator and write the gap equation down. The results in this 
section are the basis for the perturbation theory, which treats subdominant 
terms ignored at the mean-field level. 



4^1 Preliminaries 

In (l+l)-dimensional models, the chirality (i.e. the eigenvalue of 75) character- 
izes the moving directions of particles. In higher dimensions, the corresponding 
7-matrix is not 75, but r^5 = 7o7i|| for particles moving to the x^y -direction. 
We define 

1 + r 

^i± = — ^^i • (33) 

This r^5 satisfies the following algebraic relations: 

(r,5)' = i, {r,5,%o} = o, {r,5,7iii} = o, [r,5,7^±] = o. (34) 

The free Lagrangian can be decomposed into two pieces. The longitudinal part 
is defined with (+,+) or (— , — ) combinations; 

£51° = 4+ Kdo - dm)iJi+ + 4_ 1(^0 + di\\)i^i- , (35) 

and the part made of (+, — ) combinations is 

£^1" = iJi+ i ^^iJi. + iJi. i ^x^i+ . (36) 

Below, we will drop off the index i as far as no confusion arises. 

At finite density with the Fermi momentum Q, it is natural to measure mo- 
menta of fermions from the Fermi surface. Accordingly, we take fields with 
shifted momenta. 



^{x) = e'^^^w'^^^x) = e^^^^ii7/;^(x) (37) 



or 



i^±{Sp) = MSPW ± Q.P±) , {Sp= (5p||,p^)) , (38) 

in momentum space. We use the notation Sp to emphasize that momenta of 
ijj^ field are measured from the Fermi surface. Using the ip^ field, one can easily 
identify dominant and subdominant terms at large density. 
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In the T/^'-representation, the longitudinal part becomes ^^ 



4" ^ ^+ [^(^0 - 5||) + Qj^V + C [i(5o + 5||) + QJC , (39) 

and the transverse kinetic term and the mass term acquire the oscillating 
factors, 

C^l"" -^ ^'^ i ^^^'_ e-2 ^^^11 + ^'_ i ^^^'^ e^ ^^^11 . (40) 

Such oscillatory terms are suppressed near the Fermi surface by powers of 
l/Q. In the free theory, in fact, the excitation energy at \Sp\\\ <^ Q is 

e'-'^Sp) = vW+W+S -Q = \SP\\ I + ^^^ + ■ ■ ■ • (41) 

Terms with oscillating factors define what we call "subdominant" terms. 

The four-Fermi interactions can be also separated into dominant terms and 
subdominant terms, 

im' = l{{m' + (^iTsV^)') + ^((^V^)' - (V^iTs^)') 
= 2 (^+^_)(^_^+) + ((^-^+)^ + (^+^-)') 
-^ 2 (7/5;V^^)(V5'_7/;;) + (7/5'_V;;)2e^^^^ii + (7/5;7/;^_)2e-4^^"ii . (42) 

The first term corresponds to the continuous symmetric part, which becomes 
IR dominant at high density. We will apply the mean-field Ansatz to dom- 
inant terms, while subdominant terms are treated as perturbation. In this 
treatment, a gap will be found only near the Fermi surface. The gaps will 
not open periodically in momentum space because of the absence of diflFerent 
harmonica^n 

Finally, for later convenience, we write Eq. (42) in momentum space including 
the form factor explicitly. For (H — )(— +) combinations of the dominant part. 



^^ In the grand canonical ensemble the basis '0' with Q — fiq eliminates the chemical 
potential term reflecting that in the '^'-representation we can deal with dynamics 
near the Fermi surface as in vacuum. This simple logic is not directly applicable in 
the canonical ensemble since the density constraint does not explicitly appear at 
the Lagrangian level. 

^^ The current problem is different from the problem of the Peierls instability with 
an external periodic potential. In our case the coupling between the mean fleld and 
particles depend on the zb combination. For instance, the mean fleld (t/j-^t/j-) = 
^g-2iQx|| ^^^ scatter the particles from the + region to the — region, but cannot 
scatter from the — region to the + region. In this way particles and holes are kept 
around the Fermi surface, without going to the higher harmonic regions. 
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we have 

= I U\{Sp + q- 2Qn) iP'_{5p)) U'_{Sk) iP\{5k + q- 2Qn)) esp,sk • 

Jq,6p,6k ^ ^ ^ ^ 

(43) 

Note that ifqc^ 2Qn^ aU fields can be close to the Fermi surface simultaneously 
when 6p ^ Sk ^ 0. This is the reason why the configuration with q ^ 2Qn 
becomes dominant in the path integral. For this reason, we should choose the 
wave vector of the chiral spirals to be g = 2Qn in the high-density limit. This 
expression also indicates that the exciton-type condensation (i.e. homogeneous 
chiral condensation) with g = is not favored energetically. 

On the other hand, for (H — )(H — ) combinations of the subdominant part, we 
have 



/ (t/5+(p + q)iJ.{p)) {iJ+m-ik + q)) Op,k 

J q,p,k 

;(5j9 + q- 2Qn)^'_{5p)) (ij'^{5k)^'_{5k + q + 2Qi 
Q[k\ - {5p- 5k - 2Qnf) . (44) 



In contrast to the dominant terms, subdominant interactions require that at 
least one fermion must go far away from the Fermi surface regardless of any q. 
The propagation of such a fermion serves the 1/Q suppression in the quantum 
corrections. 



4^2 Formal Treatment: Bosonization 



Now let us introduce the mean-field Ansatz for one patch. The auxiliary-field 
method may be applied as before. There are two slight modifications on the 
standard approach. One is that our boson field is introduced as a complex 
field since we need to eliminate quark bilinears, ip^ip- and ip-ip^^ which have 
complex phase factors that are opposite to each other. Another is that boson 
fields are used to be replaced with only dominant four- Fermi interactions. 

Inserting an identity to the original partition function, we introduce the bosonic 
terms, 






fdx' [ ^\q',p)e,^,^q',k), (45) 

to the original action. We replace the dominant four-Fermi interaction with a 
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Yukawa- type vertex by shifting the boson iields, 

G 



^ G - ^ ^ 



$(g;fc) - 
then the Yukawa vertex is 

5^,^ = - / da;° / [ / Op^k $(g; k) ] Mp) ^-iP + O) 



- fdx' [ \ fep,k^\q;k) 

J ^p,q L ^k 



(46) 



(47) 



where terms inside of [• • • ] in the first and second terms are the self-energy, 
^ln{-Q'^P) and T.rn{q]P), respectively. 

4>3 Mean Field for Chiral Spirals and the Quasi- Particle Spectrum 



According to the arguments around Eq. (43), in the high-density limit the 
Ansatz for chiral spirals in the i-th patch is 

$o(g; k) = S{q- 2Qn,) A{k) , $^(g; k) = S{q- 2Qfi,) A{k) , (48) 

where A is a real field which characterizes the magnitude of the condensate. 
We will give what would happen if we chose different q in Appendix A. The 
Fourier transformation with respect to the total momentum of the boson fields 
gives the expression. 



^o{x]k) = A{k) e 



210x11 



^l{x]k) = A{k) e 



-210x11 



(49) 



from which we can see that the complex nature of the boson fields is taken 
into account in the phase factor. 



For later convenience, let us define the mass gap function. 



M{p) = / e,^, A{k) , 

Jk 



(50) 



which will be determined self-consistently after constructing the mean-field 
propagator. Using this shorthand notation and shifting momentum p — )• p — 
Qfii, the mass vertex becomes 



«5$,^ — )• S^ 



--- j Ax^ j [m{p- Qfii) i)+{p-Qni)i).{p + Qfii) 
+ M{p- Qfii) xp-{p + Qfii)xp+{p- Qfii) ) 
dx^ I (^M'{Sp) i)'^{5p)i)'_{5p) + M'{5p) ^'_{5p)^'^{5p) ) . (51) 
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In the last line of the above equation, we have replaced the loop momentum 
^with 5p and have defined M'{5p) = M{p — Qfii). 

Here it is very important to notice that although our inhomogeneous conden- 
sate breaks the translational invariance, momenta measured from the Fermi 
surface, 5p^ is a conserved quantity. This is true as far as we include only 
dominant terms. The conservation of 5p is violated by corrections such as 
transverse kinetic terms or subdominant parts of four- Fermi interactions, but 
they are subleading effects suppressed by powers of 1/Q. 

At leading order, thanks to the conservation of 5p^ the eigenvalue problem 
is diagonal in momentum space. Then we can formally derive the mean-field 
spectrum of quasi-particles using the mass function M\5p). The eigenvalue 
problem for the longitudinal plus mass terms is 

where a notation, ^!\5p\\) = \^'j^{5p)^ilj^_{5p)j , is introduced. The eigenvalue 
has upper and lower branches, 

EMYm = Q±^m. ( oj{8p) = ^5pl + M'{8p)^ ) , (53) 

where the gap opens at \p\\\ = Q, and its infiuence survives up to distance 
^ Ac from the Fermi surface. (For the case that the chiral spiral has a wave 
vector Q' ^ Q^ see Appendix A.) Since the energy level for particles with a 
mass gap is pushed down as compared to free particles, the energy gain inside 
of one patch is ^ M x (phase space) = M x AcQtanB. Also, note that the 
phase space does not change before and after the formation of the mass gap, 
so that the Fermi volume conservation is automatically satisfied. 

It is important to specify the relation between these upper and lower branches 
and the respective momentum regions. They can be summarized as 

E^ = Q^uj{bp\\) (P||>g), E^ = Q-uj{bp{) (P||<g), (54) 

for states moving to the + direction, and 

£;^ = g + cj(5pii) (p||<-Q), E\^^Q-^[bp\\) (p||>-g), (55) 

for states moving to the — direction. See Fig. 11 for a graphical summary. 

In Appendix A, we repeat calculations for q ^ 2Qn^ and explain that at 
sufficiently high density, the choice q = 2Qn is the best way to minimize the 
single-particle contributions with the fixed particle number constraint. 



24 






\ ^^1 

sek^\ /se^ 


dp 




^Ey/ 







Fig. 11. The mean- field single-particle dispersion in the presence of the chiral spirals. 
The particle orbit is occupied up to |p||| < Q. The Fermi volume conservation is 
satisfied before and after the formation of chiral spirals. (Left) The energy £'mf in 
the '^-representation. The gap opens at the edge of the Fermi surface. (Right) The 
energy 6E = ^'mf — Q in the '^'-representation. 

4^4 The Mean- Field Propagator in the Canonical Ensemble 



Let us construct the mean-field propagator for the quasi-particle. We are work- 
ing in the canonical ensemble with a fixed shape of the Fermi sea, i.e. one 
wedge, filled with quasi-particles up to \p\\\ < Q. This can be restated as the 
condition that all occupied states must take the lower energy branch of the 
mean- field spectra, E = Q — uj. To take into account such information in the 
propagator, we have only to take a proper ie prescription as we usually do in 
the vacuum case. 



The relevant fermion bilinear part of the action (longitudinal plus mass terms) 



IS 



ll+M. 



dpo f ^\po. Sp) ( {po - Q)i' + 5p||7" - M\5p) ) *Xj9o, ^ • (56) 






Note that Q couples 7°, not to 7II , as a consequence that we measure momenta 
of il)± in an opposite way. Defining 5pQ = Po — Q, our mean-field propagator 
in the region of condensation is, 



Smf{Sp) = i 



^Po7° + Sp±l^ + M'{5p) 
{5pof - {Sp±y - M'^{5p) + ie 



0±{Sp) 



(57) 



and out of the condensation region, it is just a free quark propagator. Here 
ie is introduced as in the vacuum case in order to classify the propagations 
of particles in the upper and the lower energy branches. The last function 6± 
takes into account the phase space restriction of one wedge, 



e^{5p) = e (\Q + 5p\\\ tune - \p±\) 0{\Q-Sp\\\ta.ne-\p, 
~ ^(gtane - \p±\) ( for g > \Spii\ ) . 



(58) 
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This means that both quasi-particles and holes must be within one wedge 
to create chiral spirals. For the most part, we will assume the approximate 
expression in the second line, provided that Q is sufficiently large. 

In the following, we will frequently use the component expression which is 
defined as 



Smf{Sp) 






/l + Ts 



2 



l-r5 l + Ts ^ l + Ts \ 

MF ^— ^— <Smf ^— 1 

i-r5 



5mf 



l-Ts ^ l + Ts 

^MF 



2 / 
(59) 



Some useful component expression is given also in Appendix B. 



4^5 The Gap Equation and Its (1+1) Dimensional Character 



Since we have the formal expression of the propagator, we can now write the 
gap equation explicitly. The equation is (with the trace for the Dirac indices) , 



M\5p) 



-G 



9, 



= 2G 



rddko r 

J 27r Jsk 

dSkii d5k± 



6p,6k 



tr 



{-i)M'{Sk) 



I (Sko) 
M'{5k) 



{2ny 



SB, + M'^iSk) 



iSk±y-M'^{Sk) + ie 
-0p,,9^{5k). 



e±{Sk) 

(60) 



When we treat the self-consistent equation, we have to distinguish two situ- 
ations depending on the fermion momenta. One case is that p± resides suffi- 
ciently far from the patch boundary. The other case is that p± is so close to 
the boundary that its mass gap is affected by other patches. 

In this section, we focus only on the former case without the boundary effects 
and postpone the discussion on the latter case to the later sections, which 
actually goes beyond the one-patch problem. 

For pj_ sufficiently far from the patch boundary, the restriction 9j_{Sk) is au- 
tomatically satisfied by the other restriction 9p^k^ so the former does not play 
any essential role. Then the gap equation has (l+l)-dimensional solutions. 
Indeed, we can find the mass gap function independent of p±^ 

M',{5p) = M',{5p\\) . (61) 

When we look for such a solution, we can factorize the integral over k^, which 
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gives the (l+l)-dimensional gap equation, 

where the p^ dependence disappears from the RHS as it should be. 

Let us note that this factorization could not yield the p^-independent solutions 
if we did not ignore corrections from the transverse kinetic terms ^ p\/Q or 
0±{6k). If one of such terms become relevant, then the p± dependence of 
the RHS does not disappear in the above manipulations. This means that the 
(l+l)-dimensional mass function can be obtained only if Q is sufficiently large 
and p± is not too close to the one-patch boundary. 

Now let us show that the non-trivial gap arises from the IR effects near the 
Fermi surface. We shall consider the case with Sp\\ = 0. If we separate the 
integral region of the RHS into regions below and above a certain scale cAf such 
that the momentum dependence in Mq can be ignored below cAf. Assuming 
that cAf < Af, we have 



Mn(0) c^ / ^ , "^ ^ + (finite positive terms) 

'^ ^ 7T Jo 27r ^Sk^ + M^^O) ^ ^ 

4GAf / cAf \ 
= Mq(0) • In + (finite positive terms) , (63) 

where the logarithmic term comes from the (l+l)-dimensional character of 
the equation, and expresses the IR effects. 

In the same way as the Cooper instability in superconductivity, we can find 
the solution of the gap equation from the IR structure of the equation. Indeed, 
if we take Mq{0) too small, the logarithmic part is divergingly large and the 
RHS exceeds the LHS much. Thus Mq{0) must be taken substantially large 
until the IR contributions are tempered to be the same order of the LHS. 

We point out that one oversimplification in the above expression is related to 
our approximation to ignore the transverse kinetic terms. If we recover those 
kinetic terms, the logarithmic part must be modified effectively as 

which tempers the growth in the IR region. Thus, the gap solution might not 
be found unless the density is sufficiently large. 
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Finally, as a solution of the gap equation, Mq is parametrically given as 

M^ - Af e-^/^^f . (65) 

where C is some number. 

5 Perturbation Theory with Chiral Spiral Mean Fields 



In this section we develop systematic computation of the corrections from 
subdominant terms up to the leading order of the 1/A^c expansion. 

Using the stationary phase approximation at large N^^ the fermionic partition 
function under chiral spiral backgrounds for fixed G is 



^^ [(*).©] = jWVii)' e^(^MF+A5) 



^2 



1 

2!^""^^ ' /mf' 



.MF\l + iA<S + -(A5f + ■••>_, (66) 



where (• • •)mf is the expectation value when we use the mean-field weight, 
q^Smy ^ in the path integral. The action is made of 

i 

as = y: {s^ + <sr'''^*) + E <5.. + E ^^.>^ + E ^^^M > (67) 

where Smf is the mean-field action which is the (uncorrelated) sum of one- 
patch actions. The S^ and 5|^^-^^* are subdominant terms inside of the i-th 
patch, which were not treated in the last section. Finally 5^^, • • • describe 
the interactions among different patches. 

The energy density functional is given by ^[($), B] = — ilnZ/Va = ^mf + Af 
(with V3 being the (2+l)-dimensional space-time volume), where 

^MF = ^ In ^MF , Af: = ^ • (i A5 + ^(A5)^ + . . . );7 . (68) 

We wish to measure the energy benefit from the chiral spiral formations. To 
do so, we need to subtract the energy of the trivial configuration, and so we 
should compute. 

When we apply the perturbative expansion soon later, it is convenient to 
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reorganize the above expression into 

ssm^e] = (£[($), e] -£[0,6]) + (£[0,6] -£[0,0]) . (70) 

The first term expresses genuine condensation effects, while the second term 
comes from the deformation energy, which was already computed at the in- 
troduction and turned out to be ^ p^Q^. Our task below is the computation 
of the quantity inside of the first parentheses. 

In the following, we will use the T/^'-representation. The advantage of doing this 
is that the momentum Sp in the propagator is conserved, and we can express 
the propagator as a function of relative distance in space-time, x — y. Noting 
this fact, many of diagrams can be easily cast away. Let us recall that the 
subdominant terms have oscillating factors multiplied to the fermion fields. 
Non-zero contributions remain only if the combination of vertices has the 
oscillating factors in a form of e^^^^^~^^ii , Q^'^^Qix-y)\\ . . . because propagators 
are functions of x — y. In other words, this is simply a consequence of the 
momentum conservation in our shifted variables. 

For example, let us see the first-order expansion, that gives a tadpole contri- 
bution, such as 

d'x{^{i,'^i^^i,'\^ e-^^-^^ii + {i''Jp±i^'^)^^ e^^^^ii) . (71) 

This is, however, a space-independent quantity times an oscillation factor, and 
so its spatial integral vanishes for Q 7^ 0. This is an example of the momentum 
conservation. Non-zero correction terms start to arise from the second order 
of the perturbative expansion. 

In the following, we will first compute the corrections from one patch which 
include transverse kinetic terms and subdominant four-Fermi interactions in 
one patch. Second, we will compute interaction terms including patch-patch 
interactions. 



6 Corrections from One Patch 



We consider the perturbative corrections inside of the i-th patch. For the 
moment we will omit the subscript i. The sources of second-order corrections 
are enumerated as follows: (1) (transverse terms)^, (2) (four-Fermi interaction 
terms)^, (3) (transverse terms) x (four- Fermi interaction terms). The last one, 
(3), cannot cancel oscillating factors out, and only (1) and (2) contribute to 
the free energy. 
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5p\\ + 20 




Fig. 12. A diagram for the transverse kinetic terms in the i-th patch. Here we 
used the '^'-representation in which '0' fields effectively feel incoming and outgoing 
momenta, 2Qni. This makes at least one fermion go far away from the Fermi surface. 

6. 1 Product of Transverse Kinetic Terms 



The non-zero contributions from the product of transverse terms arise from 

(H — )( — h) or ( — h)(H — ) combinations, i.e. 



^-^^tra.ns 



1 



d^xd?y ( i)'M^ij'_{x) i,'_i^^ij'{y) 



-2\Q{x\\-y\\) 



MF 



V3 

^"J i27ry''^J 27r 

X (-i)tr S—{Spo,Spii,p±) S++{Spo,Spii +2Q,p±] 



d^5p 2 f d(^Po 



I 



(72) 



After the po integration, we have 
dSp\\ 



Ac+rans — -l^ r. 



dp± 2 



F{5p\\,2Q) 



(73) 



where the integral of the transverse momentum is restricted within one patch. 



dp 



L«,ir-/ir<'-<*"^-"'-<*'+2«-^J 



i(Q) 27r J 27r 

and we defined the function F{5p\\, 2Q) as 

6p\\ 5p\\ + 2Q 



F(5p||,2g) = l + 



uj{5p\\) uj{5p\\+2Q) 



(0 < F < 2) . 



(74) 



(75) 



Now let us analyze Eq. (73). The key point is that at least one of the momenta, 
5p\\ or 5p\\ + 2(3, must be much larger than Ac. If both of the momenta are 
much bigger than Ac, then the mass gap does not exist, so that the result is 
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just reduced to the one in the free theory. It is thus a non-trivial case when one 
of the momenta is close to the Fermi surface. Let us consider such a situation 
in the following. 



Supposing that l^py | < A^ the other momentum ^py + 2Q satisfies 2Q — Ac < 
6p\\ + 2Q < 2Q + Ac- Therefore we can use the free fermion dispersion for 
(jj{6p\\+2Q) = 6p\\ +2(5, and can apply an approximate upper bound, Q tan B, 
for the p^-integral. 



In order to make physical interpretations, we should notice that generically the 
perturbative expansion contains the trivial contribution that is independent 
of the condensate, and this must be subtracted. Here we subtract such trivial 
contributions with the same phase space as the non-trivial one. Outside of the 
condensation domain, two contributions cancel out. After subtracting terms 
at M^ = 0, the non-trivial contribution of Eq. (73) is given by 



Nr. r^^ _ /-gtanB ^ p^ 



7V r^^c r^ tan fc> 

A^trans - ^ / d5p\\ / dp_ 

TT^ Jo Jo 



2Q 



Sp\\ 



5pl + M^^{5p\\)^ 



N^ • M'Q^ tan^ 9 {M' - A;) , (76) 



where —2 in the first line represents the subtraction of the trivial contribution. 
Note that the energy correction before performing the integration is ^ p\/Q^ 
and so at the level of the computation of the fermion dispersion relation, this 
terms were negligible effects. 



The sign of Eq. (76) is negative, which means that this term corresponds to 
an energy gain. This term can be regarded as a controllable correction only if 
it is sufficiently smaller than the energy gain from the one-patch mean field, 
A^c • Af QtanB. This requires the condition, B <C (Af/Q)"^^^, as already stated 
at the introduction. 



Finally let us mention on how the perturbative analysis of the correction from 
the transverse kinetic energy can be organized in a systematic way. If we go to 
higher order of the transverse kinetic terms, we have terms of higher powers of 
p\lQ in the integrand. Computing the p^-integration with the integral region 
^ QtanB results in terms of higher powers of B <C 1, which appears as an 
expansion parameter in the systematic computations. 
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Fig. 13. Contributions from the irrelevant vertices. The ^'-representation is used 
and all fields belong to the i-th patch. (Left) The 0{N^) contribution which is 
ignored in this work. (Right) The 0{Nc) contributions. The diagram can be inter- 
preted as the condensate-condensate interaction mediated by particle-hole excita- 
tions. One of fermions in the internal loop must go far from the Fermi surface. 

6.2 Product of Four- Fermi Interaction Terms 



We will compute the second-order perturbation of irrelevant four- Fermi inter- 
action terms, which reads 






!■ 



d'xd'y((ij'ij'_{x)) U'_fJy 



-4ig(x||-y||) 



MF 



(77) 



Equation (77) yields several distinct terms depending on how fermion lines are 
contracted. Here we introduce the hierarchy for different contractions by ap- 
plying the 1/A^c expansion. At large Nc, the dominant contributions come from 
condensate-condensate interactions mediated by virtual quark-hole exchange. 
This situation is schematically described as 



a 



2 {lp+i^'-)MF 



Nr. 




(Tztr 



S++iz)S..{-z) 



-4iQz« 



{i,'_i,') 



MF 



K. 



(78) 



with z = X — y (see the right panel in Fig. 13). Here the form factor is not 
exphcitly written yet. This contribution is 0{Nc) and positive, according to the 
previous calculations which treated the integral in the bracket. One fermion 
pair at x and y is contracted at one space-time point and becomes condensates. 
The remaining part is for the particle-hole propagation between x and y. 
Loosely speaking, this part can be interpreted as the propagation of meson- 
like objects with total momentum 4(5 (though ladder- type resummation is 
necessary to construct the meson propagation in fact). 

Let us note that the size of the condensate of 0{Nc) compensates for the 
suppression factor of 0{1/Nc) in the intrinsic interaction vertex. Other con- 
tractions without the condensate cannot be accompanied by a fermion loop of 
0{Nc) and are suppressed by l/Nc. 
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We should not take Eq. (78) literally, however, since the loop integral is UV 
divergent, which is regulated with the form factor. It should be mentioned that 
such apparent UV divergence couples to the condensate, so the subtraction 
of the trivial configuration without the condensate cannot regulate the UV 
behavior. 

Including the form factor explicitly and carrying out the integral in Eq. (77), 
we have 

N^G^ f e(A^^ - {5p- Sk + 2Qnf) e(A^^ - {Sp- St+ 2Qnf 

JSp,Sk,Sl ^ ^ ^ 






(79) 



or equivalently, 

N, I M'(5p + 2Qn) F{5p\\AQ) M'(5p + 2Qn) . (80) 

J5p ^ > cj(5p||)+cj(5p|| +4Q) V "^ ^ ) ^ ^ 

Here the phase-space restriction for the transverse momentum is not explicitly 
written. The condensate takes a finite value only if l^fcy |, \^l\\ | < Ac. Once they 
are restricted within such a domain, bp is also restricted around —2Qn^ and 
so the integral over the phase space has a UV cutoff. The integrand itself 
is suppressed by l/Q. Finally the energy cost from this contribution can be 
estimated as 

Ai 
(one-patch energy cost) ^ N^C^ -^ • AfQtanB ^ A^cAf tanB . (81) 

We see that the above has a correct dimension with the mass dimension of 
G given by negative one, Af^. Here, a phase-space factor ^ AfQtanO arises 
because one spatial momentum is not restricted by the 9 function constraint 
in Eq. (79). As expected, this contribution is much smaller than the one-patch 
mean-field contribution ~ A^c • A^QtanG, so can be safely ignored at large 
density. 



7 Contributions from Patch Boundaries: Inter-Patch Effects 



So far, we have ignored the interactions among fermions belonging to different 
patches. Due to the form- fact or effects, such interactions occur only near the 
boundaries of patches. We will discuss the impact of such effects. 

The outline for our computational procedure is the following. The conden- 
sation effects generate the following three types of the energy contributions. 
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Fig. 14. The domain where the chiral spiral mean field in one patch couples to 

V 



particles (holes) in its nearest neighbor patches. The size of such a domain is ~ A? 



due to form factor effects. 



4o„d. ~ N, (4ir'(^o, e) + £:^-p^*^^Mb, e) + ff^f '"(^b, e)) , (82) 

where c^^^^^^ and t^ are one-patch condensation energy m the region tar 

from and close to the boundaries, and ^^t*^ ^^ represents the patch-patch in- 
teraction at the boundaries. The subscript B means the boundary, and Mq and 
Mb can be considerably different. The energy density schematically behaves 
as 

^cond./A^c - ^ (-MoAf(gtane-Af)-MBA^ + .SPf^^^(MB,e)) 

~ -MoAfQ + ^ (Af2(Mo-MB)+£:Pf'"(MB,e)). (83) 

We will show that ^^f ^ ^^ is positive (at least at the level of the second- 
order perturbation theory). Its strength is determined by the size of Mb, and 
vanishes when Mb -^ 0. 

Equation (83) can be understood in twofold ways. If we regard the incoherent 
sum of the one-patch actions as our unperturbed action, Mb at the unper- 
turbed level is ^ Mq, and ^^t*^ ^^ provides relatively large, positive contri- 
butions. This means that the inter-patch interactions between different chiral 
spirals reduce the energy gain from creating the condensates at the patch 
boundaries. Instead, if we assume Mb <C Mq, then the one-patch contribution 
is Mq — Mb ^ Mq ^ Af , while f^f ^ ^^ is negligible. In both limiting cases, the 
sign of 1/B terms is positive and of 0(Af ). 

The precise estimate of the 1/B term would be given by self-consistently solv- 
ing the gap equation of Mb with ^^t*^ ^^- This is beyond our scope in this 
paper. Instead, we will give several indicative discussions to understand inter- 
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Fig. 15. (Left) An influence of the z-th chiral spiral mean- field at the patch bound- 
aries. Only the Fermi surface close to z-th patch is shown. The i-th chiral spiral can 
scatter a particle (hole) state in the (z + l)-th patch into a hole (particle) state in the 
(i — l)-th patch. Such processes are possible only for a particle and a hole close to the 
z-th patch boundary. (Right) A diagrammatic expression of the particle-condensate 
scattering. The vector p and k must be close each other. 



patch interactions at the patch boundaries, in both perturbative and non- 
pert urbative manners. 



7. 1 Preliminaries 



As we have already seen, at large A^c, the dominant correction terms come 
from condensate-condensate interactions mediated by virtual quark-hole ex- 
citations. Therefore we will take into account only such terms. Then we need 
consider only vertices in which two indices of the patches are identical, such 
as terms in £J^J or O^^^j.. 

Actually, because of the form factor, only a particle and a hole in {i ± l)-th 
patches can directly couple to the i-th chiral-spiral mean field. Therefore we 
have only to consider the following types of the vertices (see also Fig. 15), 



/ (V^i+(P + 2Qni) ipi-ip)) ijJi-i-ik) ^i^i^^{k + 2Qfii) Op^k , 



(84) 



where we have replaced a fermion bilinear in the i-th patch with the chiral 
spiral mean field. Using the shifted momentum variables, Sp = p + Qui and 
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6k = k + Q^i-i, and the T/^'-representation, Eq. (84) can be rewritten as 

/ (^:+(<5p)C('^p))Ci,-('^^)^:+i,+('^^+Q(2n^-^i-i-^i+i)) 

J6p,6k \ / ^ / 

xe[A^^-{Sp-Sk-Qn, + Qn,_if) . (85) 



In perturbative computations, it is necessary to decompose momenta of the 
fermion field into the longitudinal and transverse components. Here let us 
briefiy summarize necessary ingredients for the computations. Note that 

^i±i = COS 2B Hi =F sin 2B fiij_ , '^i±i,± = ± sin 2B Hi + cos 2B Hij_ . (86) 

If we write 6k = Sk\\Hi-i + k±Hi-i^±^ then the momentum in ip[^i^ can be 
decomposed into the n^+i and the n^+i,± directions, 



Sk + Q{2Hi - Hi_i - Hi^i) 

Sk\\ cos 4B + k± sin 4B + (3(2 cos 2B - cos 4B - 1) 



rii+i 



+ 



6k\\ sin 4B + k± cos 4B + Q{-2 sin 2B + sin 4B) 



^i+i,j 



^ {5k\\ + 4k^e) Hi+i + fc^ n^+i,± + 0{5k\\e, QB^ 



(87) 



where we did not explicitly write quantities which are much smaller than Af. 
Similarly, let us simplify the expression of the argument in the form factor. A 
decomposition of the momentum, 

Sp — Sk — QHi + QHi-i 

5p\\ — 5k\\ cos 2B — k± sin 2B + Q{co8 2B — 1) Hi 



+ 



p± — k± cos 2B + 5k\\ sin 2Q — Q sin 2B 



^i+i,± 



- {5p\\ - Sk\\) H, + {p^ -k^- 2Qe) H,^ + 0(5fc||B, QQ^) 
leads to 



(89) 



The 6 function is non-zero only around the patch boundary, that is, p^ ^ Q© 
and k±_ ^ —QQ. 

Finally let us note that the projection operators in one patch are different 
from those in other patches. Indeed, 

1 + Ti+i^s 1 - ri_i,5 _ 1 + 7o7i+i,|| 1 - 7o(7i+i,|| cos 26 - ^i+i± sin 26) 



( (1 - cos 26) + sin 267o7i+i,± ) 



1 + r 



i+1,5 



4 



(90) 
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and 



( (1 + cos 26) T sin 26 7o7.+i,± ) • (91) 



In what follows, these slight modifications provide only negligible contributions 
of 0(6), so we need not care them seriously in the computations. 

1.2 Perturbative Consideration 



The purpose of this subsection is to get the typical size of patch-patch interac- 
tions, and more importantly, to investigate whether interactions are attractive 
or repulsive. The latter can be done even without detailed estimates of the 
gap near the patch boundaries. 

The perturbative computations proceed in almost exactly the same way as 
before. After taking the residue, we have an expression analogous to Eq. (80), 

N, f M-{5k,M + 2QQ) ,,,^"^^^";tt^^L n^ > 0' (92) 

Ait ^^ " ^ ujB{Sk\\)+ujB{Sk\\+Akx_Q) " ' ^ ^ 

where < Fb < 2 by definition of Eq. (80), and k^ ^ —QQ- The product of 
Mg comes from the condensates (i.e. fermion loops) in the i-th patch, and the 
remaining piece comes from virtual particle-hole excitations. We emphasize 
that for the momentum conservation to be satisfied, both condensates must 
come from the same patch. The subscript B is attached to remind that the 
sizes of the gap and mass function near the boundaries may be considerably 
different from those far from the boundaries. 

The sign is positive, so this is an energy cost. Although this is nothing beyond 
a generic fact inherent to the second-order perturbation, it indicates that 
boundary effects tend to temper the magnitude of the condensates near the 
boundaries. 

It is very important to notice that the energy cost may be comparable to the 
energy gain from forming a condensate within the same phase space ^ Af . In 
contrast to the case of subdominant terms, we do not have 1/Q suppression 
because all fermions can move around the Fermi surface during the virtual 
processes. 

Let us investigate the order of the magnitude. If we assumed Mg ^ Q6^, we 
can make an approximation that 

Integrand in Eq. 92 ^ ^^ . \ . 1 - ^^ " ^^^,, + • • • , 93 
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thus the integral of the above integrand in the IR region gives an approximate 
expression as 



AfM^2 



/Af\ 2ge2 

m 



MJ MTm§ 



- M^ (^ + o(ge^)j, (94) 



where we have used the parametric behavior of the mass gap, Mg ^ Af e~^/^^f . 
This expression indeed confirms that the energy cost is of the same order as 
the energy gain from condensation effects within the same phase space. 

Unfortunately, any obvious expansion parameter did not appear for this per- 
turbative expansion, so we have no good reason to cast away higher-order 
diagrams. For a more reasonable estimate, the non-perturbative computa- 
tions are necessary to simultaneously treat patch-patch interactions and the 
mean- field problem near the patch boundaries. 



7.3 Some Non-perturbative Considerations: A (1 -hi) -Dimensional Example 
of Two Chiral Spirals 



To get some insights for the inter-patch interactions between several chiral 
spirals, let us consider the simplest (l+l)-dimensional example. We assume 
the mean field which has two chiral spirals with wavevectors Qo and Qi. We 
will see how these two chiral spirals affect each other. 

The mean-field Hamiltonian is 

^MF Hz) = [^^ ^_.Q,. ^ ^^ ^_.Q,. .^^ j Hz) , (95) 

where ^{z) = f'0+(z), '0_(z) j defined as before. 

We rewrite this expression in the '^'-representation in order to eliminate the 
oscillating factors. When we have two chiral spirals, those oscillating factors 
cannot be eliminated simultaneously. If Mq > Mi, it is better to eliminate 
^±2iQoz^ as will be clear in the following. Using t/;^ = ipe^^^^ and multiplying 
the Hamiltonian squared, we have the following Schrodinger equation, 

where 5Q = Qo — Qi, and 

^diag. = -dl + (Mo - Mi)2 + 4MoMi cos^ 5Qz . (97) 
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Fig. 16. The potential V{z) — Q. (Left) SQ = case. The potential is constant. 
(Right) 6Q y^ case. The potential is oscillating, so particles can stay around the 
valley of the potential. For very small 6Q, its kinetic energy cost is small. 

The off-diagonal Hamiltonian has the amplitude proportional to SQMi. There- 
fore, if SQ or Ml is small enough, one can ignore the off-diagonal part. Here 
one can understand that our choice to eliminate Qq rather than Qi is suited 
for this approximation since Mq > Mi. 

The diagonal part has a positive oscillating potential whose period is l/SQ] 
see Fig. 16. A singularity lies dit 5Q = which leads to a constant potential, 
and the energy spectra are discontinuous from those at 5(5 7^ 0. If SQ is small 
enough but not zero, then we can find the eigenfunction whose kinetic energy 
is very small. That is. 



£;- Q± |Mo-Mi|, 



(98) 



for < |5(5| <C Mq. When Mq is comparable to Mi, they nearly cancel each 
other, making the effective gap small. This analysis indicates that two chiral 
spirals with different but similar wave vectors tend to reduce the energy gain 
in single particle contributions. 

If two chiral spirals have substantially different wave vectors, say, SQ ^ M^ 
then the valley of the potential becomes narrow, and the kinetic energy is 
^ M. In such a case inter-patch interactions do not strongly reduce the energy 
gain. This remark will be increasingly important when we consider the lower 
density region where higher harmonics of chiral spirals start to contribute 
because of the subdominant terms. We will discuss this in the next section. 



8 Discussion 



So far we have used the four- Fermi interaction which is strong enough to induce 
the chiral symmetry breaking near the Fermi surface. At the same time, we 
have relied on the high-density approximation to make our discussions simple 
enough to explore analytic insights. In reality with A^c = 3, however, gluons will 
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be screened at high density, reducing the magnitude of our couphng constant. 
Perhaps our approximations may need improvements to be reahstic. 

Hence we try to extrapolate our insights at high density into the lower-density 
domains, by arguing how correction terms grow up. Essentially interweaving 
chiral spirals are disturbed by the transverse dynamics and inter-patch interac- 
tions, and these effects become increasingly important at lower density. In this 
respect, we quote several elaborated numerical studies in the low-density side, 
in order to complement our current studies. At the same time, we will make 
use of our perturbative corrections to interpret some results in the existing 
literatures. 

Another important question is how the 1/A^c corrections grow up at higher 
density. We will summarize effects which we have ignored in the large-A^c 
approximation. 

The remaining topic is the instanton-induced interaction [27], which is fre- 
quently used to introduce strong diquark correlations [28] . We will see that this 
interaction would provide quite different effects below and above the strange- 
quark threshold. 



8.1 Comparisons with Other Works in the Low Density Regime 



In Ref. [29], the authors have numerically studied the chiral crystals in (3+1) 
dimensions in the density region, /Xq = 0.4—0.6 GeV, using the NJL model and 
the model with the instanton-induced interactions. They studied the scalar- 
isoscalar channel with the plane- wave oscillations ^ ae^^'^. In the NJL model, 
the strength of the four- Fermi interaction becomes weaker at higher density be- 
cause of the model cutoff, so they found only small mass gaps, ~ 0(10) MeV. 
On the other hand, the instanton-induced interaction is stronger near the 
Fermi surface, so it is possible to have a large mass gap of O(Aqcd)- Impor- 
tantly, they showed that the creation of the differently oriented chiral-density 
waves (crystals) does not provide much energy benefit, and the chiral density 
wave evolved only in one direction is energetically favored. This result can be 
interpreted as a consequence of inter-patch interactions, and is consistent with 
our current analysis. 

In the case of the chiral density wave in one particular direction, there can 
be a better solution than the simple plane wave. Recently, such a solution in 
the (3+l)-dimensional NJL model was found by the authors of Ref. [15] who 
have used trial functional motivated from (l+l)-dimensional studies [30]. 
The following results are deeply relevant to ours: (i) The solution appears to 
be of a solitonic type at low density, and approaches the plane- wave type at 
high density, (ii) A quark number modulation also occurs at low density. It 
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smoothly approaches the uniform distribution as the density increases, (iii) 
The chiral spirals with a and vf, which is naively expected, is not energetically 
favored as compared to single modulation of a, as far as we keep /lLu = f^d- 

These statements are all consistent with our analyses, and in fact could be 
inferred from our framework as follows. Let us explain this in order. 

(i) The deviation from the plane-wave solution is caused by the subdominant 
terms in our formulation, 

7/5;a/i7/;'_e-2i^^ii , (7/5;7/;^)2e-^^^^ii , • • • (99) 

which provide higher harmonics necessary to construct solitonic solutions at 
low densityp^ As we discussed, these terms become unimportant as the den- 
sity increases, recovering the plane- wave solutions. 

(ii) If we compute the expectation value of the quark number, the non-perturbative 
mean-field propagator gives uniform distribution, while the perturbations from 
subdominant terms can generate the spatial modulation. The distribution ap- 
proaches the uniform one as the density increaseq^ 

(iii) It is quite straightforward to extend our one-flavor studies to the multi- 
flavor ones in terms of $ = (t^, rf, • • • )^, and we can easily infer that the chiral 
spirals should emerge as a rotation in the U{1) quark number sector, 

($+$_) = A e-^^^^ii , ($-$+) = A e^^^^ii , (100) 

which are equivalent with the following combination, 

($$) = 2Acos2Qx|| , ($i7o7||$) = 2Asin2Qx|| . (101) 

The expectation value of the latter was not calculated in Ref. [15], though. 
Here we have not found any particular mechanism to generate flavor rotations, 
at least in the high-density limit. Of course, once we had explicit flavor break- 



^^ In (1+1) dimensions at T = 0, one can validate this discussion by investigating 
the Gross-Neveu (GN) model [31] with or without the continuous chiral symmetry. 
The former is free from the subdominant terms, and the chiral spirals can appear at 
arbitrary low density and the quark density is always uniform. (At nonzero T, we 
have the twisted kink crystal in which the amplitude field also modulates [32].) On 
the other hand, in the version with discrete chiral symmetry, solitonic objects first 
appear at density beyond some critical chemical potential which is slightly lower 
than the constituent quark mass. As density increases, subdominant terms stop to 
disturb the chiral rotation, then quark distributions smoothly approach those with 
the chiral spirals and uniform quark density. See also Sec. 6 in Ref. [6]. 
^^ Explicit calculations will be reported elsewhere. 
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ing coming from conditions such as the charge neutrahty and ;5-equihbriump^ 
they are hkely to generate other chiral rotations as weU. It would be interesting 
to consider astrophysical consequences of such inhomogeneous distributions. 
(See discussions on the imphcation to the ghtch problem in Ref. [34], for ex- 
ample.) 

Assembling these works and insights in this paper, let us infer what kind of 
calculations are desirable at low density. While the chiral crystals were not 
favored in the plane- wave Ansatz in Ref. [29] , we know that higher harmonics 
become increasingly important at lower density, as shown in Ref. [15] and our 
perturbative calculations. They are relevant to describe a localized quark num- 
ber density as well. An interesting question is whether crystals including higher 
harmonics are energetically more favored as compared to one-dimensional soli- 
tonic configuration in Ref. [15]. As discussed at the end of the previous section, 
inter-patch interactions are strong for the chiral spirals with close wavevectors. 
That argument, however, also implies that two chiral spirals with very diflFer- 
ent wavevectors do not strongly destroy one another. Thus, for configurations 
with higher harmonics, deconst ruction due to inter-patch interactions might 
be tempered, so that the solitonic crystal structure might be energetically 
favored. 



8.2 On the l/N^ Corrections 



As the density increases, the 1/A^c corrections grow up because the increas- 
ing phase space around the Fermi surface enhances low-energy quark fiuctu- 



atione ^^ . The enhancement is parametrically ^ {l^q/ ^qc'dY~^ where d rep- 
resents the number of spatial dimensions. Such effects are illustrated in the 
typical diagrams in Fig. 17. The diagrams (a) and (b) reduce the eflFective size 
of our coupling constant as the density increases. The diagram (c) modifies 
our resummation scheme or the mean-field treatment shown in Fig. 4. In the 
terminology of the four- Fermi interaction, we have treated the Hartree term 
while ignoring the Fock term in the large- A^c limit, which will be modified. 

Especially let us note that in the diagram (c), all momenta of the loop, the 
incoming quark, and the outgoing quark need not to be close, in contrast to 
the leading- Ac contributions. Therefore once the 1/Ac contributions to the 
condensate become comparable to the leading order, they will violate the 



^^ These conditions are driving mechanism to destabilize a homogeneous color- 
superconducting phases [33] into the crystalline states [34,35,36,37]. 
^^ Needless to say, when we try to include the 1/Ac corrections, we have to restart 
all of the computations from the vacuum problem, in order to renormalize the theory 
including fermion loops. Here we are arguing the medium-induced modification after 
the correct renormalization is made. 
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{a) (b) (c) 

Fig. 17. Some diagrammatic examples of the 1/Nc corrections for: (a) the gluon 
propagator, (b) the quark-gluon vertex, (c) our resummation scheme shown in the 
left panel in Fig. 4. Modifications by three- and four-gluon vertices should exist as 
well. 



locality of the quark-condensate interactions in momentum space. It means 
that inter-patch interactions among chiral spirals occur not only near the patch 
boundaries but everywhere near the Fermi surface. The situation becomes 
much more complicated than that we have treated in this paper. 

The quantitative estimation of the 1/A^c corrections is a highly non-linear prob- 
lem. It should be quite sensitive to the effective masses of quarks, while the 
mass is determined by non-perturbative interactions whose strength is reduced 
by the screening effects by quarks. In this respect, it does matter whether we 
took into account possibilities of the inhomogeneous condensates or not. Un- 
der the assumption of homogeneous chiral symmetry breaking, quarks lose 
the mass gaps originated from quark-antiquark condensation, and then fluc- 
tuations would rapidly grow up at finite density. However, if we included a 
possibility of chiral symmetry breaking near the Fermi surface, the quarks 
acquire the mass gap through the constituent quark mass, suppressing the 
fluctuations near the Fermi surface. Calculations including the latter have not 
been done yet. 

The above issue is also related to the validity of the stationary phase approx- 
imation. Here we have to distinguish two kinds of the fluctuations. One is 
the amplitude fluctuation, largely related to the size of the quark mass gap. 
Another is the phase fluctuation, which is massless. 

For the latter, one might suspect that, according to Mermin-Wagner-Coleman's 
arguments on lower dimensions [38] , the IR fluctuations of the phases would 
be strong enough to destroy chiral spirals with quasi (l+l)-dimensional struc- 
ture, and therefore quarks would lose their mass gaps. 

A few remarks are in order. First, the absence of the spontaneous symmetry 
breaking in real (1+1) dimensions is extremely sensitive to the deep IR region. 
Our phase fluctuations cannot access such a region since the Fermi surface 
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always has finite curvature effects which provide the IR cutoflp^ ^ p\/P¥' So 
even without using the large- A^c limit, we conclude that the chiral spirals have 
the long-range order. This issue was already addressed in Ref. [19]. 

Second, even if the system has only quasi long-range order, the mass gap is not 
washed out. The absence of the condensate does not mean the presence of a 
gapless quark, as far as the amplitude field takes a finite expectation value. For 
a (l+l)-dimensional example, see [39]. So a primary question concerned with 
the stationary phase approximation is whether the fields strongly fiuctuate or 
not. 

Besides the fiuct nation effects, we should also mention on the possibility of the 
diquark condensate at A^c = 3, which could not be addressed in the large- A^c 
limit. The Meissner effect would change our non-pert urbative forces. This issue 
is beyond our scope, and can be addressed only by taking the energy compe- 
tition between the interweaving chiral spirals and the color-superconducting 
phases. 



8.3 On the Six-Fermi Interaction 



Let us briefiy mention on the effects of the six- Fermi interaction that breaks 
U{l)x symmetry for three fiavors. When we discuss dynamics oi u^d quarks, 
we typically take the expectation value of ss^ and normalize the coupling 
constant as 

- G^ {uu){dd){ss) — > {G^{ss)) {uu){dd) = Gl {uu){dd) , (102) 

where we did not write any explicit structure of the Dirac 7 matrices. Thus at 
zero density, this interaction would renormalize the coupling constant of our 
model. 

At high density, however, the four-Fermi and the six-Fermi interactions show 
interesting differences. Let us first consider the following situation; fi c^ fi^ — 
jLLfi c^ fig ^ Aqcd, where all light fiavors form the quark Fermi sea. Using the 
^/^^-representation, the interaction looks like 



Q( (22>'_)(JVd'_)(5V5^_)e-^^^^^ii + (22>'_)(JVd^_)(5'_5V)e-2i^^"ii + • • 

(103) 



^^ One might think that this cutoff argument would be inconsistent with our (1+1)- 
dimensional reduction of the gap equations. This is not so. The point is that the 
mass gap is much less sensitive to the deep infrared region — even if we omit such 
small phase space, we can find the solution, like the NJL model in vacuum. 
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Here let us note that all of vertices must have oscillation factors, i.e., they are 
subdominant terms. These terms tend to disturb the chiral spiral formations, 
but eventually become negligible at very high density. 

On the other hand, at intermediate density such that only u and d quarks have 
the Fermi sea, the vertex (s^s^) is uniform and does not include the oscillating 
factor. Then we can find the dominant four-Fermi interactions which help the 
formation of the chiral spirals in the u^ d channels. Actually this is the situation 
discussed in the aforementioned work [29] . It is interesting to see whether any 
drastic changes happen near the strange-quark threshold. 



9 Summary 



In this paper we have argued how to construct the interweaving chiral spirals 
at high density. The patch size is determined by the balance between the 
energy gain from condensation and the energy cost from the deformed Fermi 
surface. The Fermi surface effects drove the spontaneous breaking of chiral, 
translational, and rotational invariance. 

A key ingredient of our discussions was that the interaction among quarks and 
the condensate, which is dominant at large A^c, was local in momentum space. 
Because of this property, energy contributions could be characterized by the 
limited phase space for the scattering between quarks and the condensate. 
In particular, complicated interplay between differently oriented chiral spirals 
happen only near the patch boundaries, so we could roughly estimate the 
energy cost by just seeing the phase space for patch-patch interactions. 

We have argued the (2+l)-dimensional model for which the geometric shape of 
the Fermi sea is relatively simple. In reality with higher dimensions, however, 
such simplicity is no longer the case since there are many ways to interweave 
chiral spirals. On the other hand, the principle to choose the best shape should 
be relatively simple according to the arguments in this paper. We need get the 
largest energy gain from condensation by maximizing the area of the Fermi 
surface until the kinetic energy cost becomes crucial, and at the same time, 
we have to minimize the length of the patch intersection lines to reduce inter- 
patch interactions. To find such a geometric shape would be a mathematically 
well-defined problem. 

The interweaving chiral spirals will survive in (3+1) dimensions, since the 
key point is the IR enhancement of the interaction, which is the fundamental 
property based on QCD. This situation is quite different from charge or spin 
density waves in condensed matter physics, which are energetically favored in 
one-dimensional systems, but not in higher dimensions [40]. 
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The form of the chiral Lagrangian near the Fermi surface was derived in 
Ref. [19], regarding the system near the Fermi surface as the (l+l)-dimensional 
chains with the transverse hopping. In the dispersion relation for the coUec- 
tive modes, transverse kinetic terms are suppressed by powers of A^/Q. As the 
density increases, therefore, the spectrum wiU approach the (l+l)-dimensional 
one, so that the IR fluctuations wiU be stronger. We also anticipate that the 
temperature effects strongly enhance the phase fluctuations [41]. Results on 
this issue will be reported in the future. 
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A The single particle dispersion for q^ 2Qn 



In this work we have assumed g = 2Qn in the (1+1) dimensional construction. 
Here we show what will happen if we chose a different wave vector. To do this, 
we go back to the original ip bases, and write the eigenvalue equation. 



where we omit the form factor effects on M for the sake of simplicity. To 
diagonalize it, we choose ip± = ip'l_e^^^ ^n , and eliminate the oscillating factors. 
The energy for the right moving branch is 
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Fig. A.l. The single particle dispersion E^ for Q' < Q (left) and Q' > Q (right) cases. 
Because of the particle number constraint, particles must occupy up to p = Q. The 
Q' — Q case with the energy E is also plotted as a guideline. 



E'^^ipw) = Q' ± Mpw - Q'Y + M^ 



(A.2) 



and its behavior for Q' < Q and Q' > Q i^ shown in Fig. A.l with the 
guideline of the Q' = Q case. With the particle number constraint, we have 
to put particles up to the momentum Q. For Q' < Q^ the particles occupy 
the upper branch, costing energy. For Q' > Q^ the particles occupy only the 
lower branch, but the energetic benefit is small because of the distance from 
the gap region. Therefore the case oi Q' = Q is the best way to minimize 
the total single particle contributions. If we add the corrections to the (1+1) 
dimensional approximations, it can deviate from Q' = Q^ though. 



B A component expression of the mean field propagator 



In the computation of the perturbative corrections, it is useful to use the 
component expressions which are decomposed into different pole structures. 
Explicitly, {uj = (jj{Sp)) 



S__(irt=ixl^x| 



l-5p\\/uj l + 5p\\/u 

Spo — u + ie Spo + u — ie 

l + Sp\\/uj l-Sp\\/uj 



Spo — u + ie Spo + CO — ie 



(B.l) 



and 



S+-{5p) = — t^Sm , S-+{Sp) = t^Sm 



(B.2) 
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where 

„ .^ . . M 1 1 



2u 



Spo — CO + ie 5po + u — ie 



(B.3) 
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